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ABSTRACT 

An empirical Compton up-scattering model is described which reproduces both the fractional ampli- 
tude (RMS) vs. energy and the soft time lags in the w 830 Hz QPO observed in 4U1608-52 on Mar. 
3, 1996. A combination of two coherent variations in the coronal and soft photon temperatures (with 
their relative contributions determined by enforcing energy conservation) gives rise to the QPO's energy 
dependent characteristics. All input parameters to the model, save a characteristic plasma size and the 
fraction of Comptonized photons impinging on the soft photon source, are derived from the time-averaged 
photon energy spectrum of the same observation. Fits to the fractional RMS and phase lag data for this 
kilohcrtz QPO imply that the spatial extent of the plasma is in the range from ~ 4 to 15 km. 

Subject headings: accretion, accretion disks — radiation mechanisms: thermal — stars: neutron — 
X-rays: stars 



1. INTRODUCTION 

Soon after the launch of the Rossi X-ray Timing Ex- 
plorer (RXTE), rapid (300-1300 Hz), nearly periodic vari- 
ability in the X-ray light curves of low mass X-ray bi- 
nary systems (LMXB's) was discovered (Strohmayer et 
al. 1996, van der Klis et al. 1996). These oscillations, re- 
ferred to as kilohertz QPO (quasi-periodic oscillations), 
have now been observed in over two dozen neutron star 
bearing LMXB's (see van der Klis 2000 for a review). 
They are distinguished by high frequencies and high qual- 
ity factors (Q=FWHM/frequency), and tend to be seen in 
pairs, with nearly constant frequency separation between 
the lower and higher frequency peaks (also called the lower 
and upper QPO). The high frequencies of kilohertz QPO 
are thought to tie them to phenomena taking place in the 
inner regions of accretion disks surrounding neutron stars, 
making them potentially valuable probes of strong gravity. 

A broad range of theoretical ideas have been proposed 
to explain the phenomenology of kilohertz QPO. In these 
models, one of the frequencies is generally identified as the 
Keplerian frequency of the innermost orbit of an accretion 
disk. The sonic point model (Miller et al. 1998) identifies 
the second frequency as the beat of the primary QPO with 
the spin of the neutron star. Stella & Vietri (1999) have 
proposed a general relativistic precession/ apsidal motion 
model wherein the primary frequency is the Keplerian fre- 
quency of a slightly eccentric orbit, and the secondary is 
due to the relativistic apsidal motion of this orbit. On 
the other hand, in the two oscillator model (Osherovich 
& Titarchuk 1999, applied to 4U1608-52; Titarchuk & 
Osherovich 1999), the secondary frequency is due to the 
transformation of the primary (Keplerian) frequency in the 
rotating frame of the neutron star magnetosphere. The 
strength (and weakness) of these models lie in their abil- 
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ity (or inability) to predict the variation of the frequency 
separation with the kilohertz QPO frequency and the vari- 
ation of their frequency with other low frequency QPO ob- 
served in the source (see van der Klis 2000 for a comparison 
of recent observations with model predictions). However, 
these dynamic models generally do not address how these 
dynamic oscillations affect, or couple to, the radiative pro- 
cesses which finally give rise to the oscillations seen in the 
X-ray spectra. The observed energy dependent features 
of the QPO, namely the strength of the QPO vs. photon 
energy (fractional root mean square, or RMS, vs. energy) 
and the phase lag of fluctuations at different photon en- 
ergies relative to variations at a reference energy, are ex- 
pected to depend primarily on the radiative mechanism 
and its coupling to the dynamic behavior of the system. 

While a unified model for kilohertz QPO, which self- 
consistently incorporates both the dynamical and radiative 
mechanisms, remains illusive, significant progress can be 
achieved by studying the radiative response of a system to 
an oscillating perturbation. Such an analysis would make 
definite predictions about the energy dependent features 
of the QPO, thereby placing constraints on the radiative 
mechanism, and perhaps, shed light on the nature of the 
coupling between the dynamical and radiative processes. 

The power law photon energy spectra often observed in 
LMXB's suggests that the dominant radiative mechanism 
in these systems is Compton up-scattering of soft photons 
in a hot plasma (Sunyaev & Titartchuk 1980, Podznyakov 
et al. 1983). The temporal behavior of such models has 
been studied in a variety of circumstances (e.g. Wijers 
et al. 1987; Hua, Kazanas & Titarchuk 1997; Bottcher & 
Liang 1999). For QPO in particular, both the fractional 
amplitude vs. energy and the relative time lag vs. energy 
can be predicted based upon models where a hot, uniform 
plasma is illuminated by soft black body photons, account- 
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ing for the system's response to variations in either the soft 
photons or the plasma (Wijers 1987, Lee & Miller 1998). 

In Compton up-scattering models, a photon's mean res- 
idence time in the hot plasma increases with the energy 
of the escaping photon since, on average, higher energies 
require a larger number of scatters. Thus it would seem 
that for these models, variability in high energy photons 
should be delayed relative to their lower energy counter- 
parts (i.e., a hard time lag is expected). In fact, analysis 
of simple Comptonization models, where the variability is 
caused only by oscillations in the intensity of the soft in- 
put spectrum, reveal that the time lag between the high 
energy, Eh, and low energy, Ei, photons should roughly 
scale as \og(Eh/ Ei) (Nowak & Vaughan 1996). It is there- 
fore rather surprising that soft lags (i.e., variability in soft 
photons delayed with respect to hard photons) seem to be 
the norm in kilohertz QPO (Vaughan ct al. 1998, Kaaret 
et al. 1999, Markwardt et al. 1999). This inconsistency 
has led to the suggestion that Compton down-scattering 
could be the dominant radiative process in these systems. 
However, it is difficult to reconcile the observed increase 
in the amplitude of the oscillations with energy with such 
models. Moreover, this would point to the unlikely sce- 
nario in which the radiative process for the time-averaged 
spectrum (Compton up-scattering) differs from that for 
the QPO oscillations (Compton down-scattering). 

This paper describes a Compton up-scattering model 
that gives a consistent interpretation of the time-averaged 
photon spectra as well as the energy dependence of both 
the fractional RMS amplitude and time lag for the 830 Hz 
QPO detected in 4U1608-52 on March 3, 1996 (Berger et 
al. 1996). The observed soft lag is a natural consequence 
of the model provided that the temperatures of the plasma 
and of the soft photon source oscillate coherently at the 
QPO frequency. It is further shown that coherent oscilla- 
tion of the temperatures is expected when energy transfer 
between the plasma and the soft photon source is taken 
into account. The model imposes limits on the size of the 
Comptonizing region and places constraints on any future 
model invoked to explain how the dynamical origin of the 
QPOs couples with the Comptonizing medium in LMXB's. 

2. THE COMPTON UP-SCATTERING MODEL 

Consider a uniform hot plasma characterized by its tem- 
perature (T c ), size (I) and optical depth r = n e axl, where 
n e is the electron number density and <tt is the Thomson 
cross section. The hot plasma up-scatters soft photons 
from an external source. Unlike previous work, the soft 
photons are assumed to be described by a Wien spectrum 
(at a temperature T s ) instead of a black body. Thus the 
soft photon source lies in a Comptonizing region with an 
optical depth large enough that the emergent spectrum 
saturates to a Wien peak. The choice of a Wien spectrum 
rather than a black body spectrum is motivated by the 
expectation that (for T s s=s 1 keV and for reasonable elec- 
tron densities << 10 25 cm -3 ), the photon density inside 
the source will not be in thermal equilibrium. It should 
be noted that broad band observations of the black hole 
system Cygnus X-l reveal that the soft photon compo- 
nent (which is assumed to provide the seed photons for 
Comptonization) is better described as a Wien peak rather 
than a black body or a sum of black bodies (Di Salvo 



et al. 2000). Since the soft photon source could itself be 
heated by the Comptonized photons from the hot plasma, 
it is assumed that a fraction, /, of the Comptonized pho- 
tons is incident on the soft photon source. In the frame- 
work of this model, oscillations in the photon escape rate 
from the hot plasma are due to corresponding coherent 
oscillations in the plasma temperature (T c ) and the soft 
photon source temperature (T s ) while the size (I), the elec- 
tron density (n e ) and the total rate of photons emerging 
from the plasma remain constant. 

For non-relativistic temperatures (kT c « m e c 2 ) and 
low photon energies (E « m c c 2 ), the evolution of the 
photon density inside the hot plasma is governed by the 
simplified Kompaneet's equation (with the induced scat- 
tering term neglected) (Kompancets 1957) 
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Here n 7 is the photon density inside the plasma, E is 
the energy, and t c = l/cr is the Thomson collision time 
scale. The number of photons per second per unit vol- 
ume injected into the plasma from the soft photon source 
is n s = CT~ 3 E 2 exp(—E/kT s ) where C is a normaliza- 
tion constant. This is balanced by the corresponding rate 
of photons escaping from the plasma n esc w n 7 / '(N esc t c ) 
where N esc w t(t + 1) is the average number of scatterings 
a photon undergoes before escaping. 

The time-averaged spectrum arising from the plasma 
can be estimated using h esc0 ~ n^/ (N esc t c ), where n 7 o is 
the steady state photon number density corresponding to 
time- averaged values of plasma temperature (T c q) and soft 
photon Wien temperature (T s0 ). For kT sa « E « kT c0 , 
it can be shown from equation (1) that the spectrum is a 
power-law with photon index (r) related to N esc by 
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Lee & Miller (1998) studied the temporal behavior of a 
Comptonizing medium by computing the linear response 
of the photon spectrum to a sinusoidally varying physical 
parameter at a given angular frequency to. Following their 
formulation, the plasma temperature is approximated as 
T c = T c o(l + <5T c e~ 4Wt ) and the soft photon temperature as 
T s = T s0 (l + 5T s e- ibJt ) where ST C and ST S are, in general, 
small complex quantities. The resulting time varying pho- 
ton number density is taken to be n 7 = n 7 o(l + 5n 7 e~ tuJt ) 
with the output spectrum given by n esc — n 7 /(N esc t c ). 
Sriy is related to ST C and ST S by the linearized Kompa- 
ncets equation (Lee & Miller 1998): 



iu>t c )n 7 o<5n 7 



-^[4£fcT c0 n 7 o((5T c + «5n 7 ) 



m e c 2 dE 
—E 2 njo5n. 
d_ 

dE 



-—E 2 kT cQ n 7Q (ST c + Sn 7 )] 
E 



-t c h s o{ 



kT s0 



3)5T S . 



(3) 



The above equation implicitly assumes that the soft pho- 
ton source spectrum h s is a Wien peak with a constant 
photon number flux. It is useful to define 8n lc as the pho- 
ton variation <5n 7 corresponding to a unit variation in the 
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plasma temperature (5T C = 1) and similarly Sn ls corre- 
sponding to a unit variation in the soft photon tempera- 
ture (8T S = 1). Thus, 5n JC (Sn 7S ) can be computed from 
equation (3) by setting ST C = 1 (5T S = 1) and ST S = 
(ST C = 0). Since only the linear response is being consid- 
ered, the variation in total photon number density can be 
written as 

<5n 7 = ST c Sn 7C + ST s Sn 7S . (4) 

Extending this method to include higher order terms leads 
to harmonics of the fundamental represented by the linear 
response. In the absence of any observed kilohertz QPO 
harmonics and in view of the simple assumptions made in 
this model, only the linear response will be considered. 

Variations in the plasma temperature will induce a vari- 
ation in the plasma luminosity. Since a fraction (/) of the 
Comptonized photons re-enter and heat up the soft pho- 
ton source, a variation in the plasma luminosity causes a 
variation in the soft photon source temperature. Thus 5T C 
and 5T S are related by 

where A, the Compton amplification factor, is defined as 
the ratio of the plasma to soft photon source luminosities. 
Using equation (4), the above equation can be written as 

fA J n l0 5n lc EdE 



ST„ » 6T C 



(6) 
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For a given variation in plasma temperature (ST C = 
\ST C \) and time-averaged model parameters (T c , T s , N esc , 
I, A and /), the fluctuations in the photon number den- 
sity (Sn-y) can be obtained from equations (3), (4) and 
(6). The complex quantity 8n 1 {E) leads to the time lag 
vs. energy and amplitude vs. energy for the QPO. Since 
the disk parameters (T c , T s , N esc and A) can, in principle, 
be constrained by fitting the time-averaged spectrum of 
the source, <5n 7 (E) depends only on the amplitude of the 
plasma temperature variation, the plasma size, and the 
fraction of photons incident on the soft photon source. 

3. RESULTS 

The predictions of the Comptonization model described 
in §2 are compared with the RXTE (PCA) observations 
of the LMXB 4U1608-52 on March 3, 1996 (Berger et al. 
1996), where a 830 Hz QPO (the lower peak - see Mendez 
et al. 1997) was discovered. The time-averaged spectrum 
from this observation is well described by the model with a 
reduced x 2 = 1-03 in the energy range 3-20 keV. The soft 
photon Wien temperature and the photon index are well 
constrained to be kT s = 1.11±0.04 keV and T = 3.55±0.05 
respectively. However, since no cutoff in the spectrum was 
observed, spectral fitting imposed only a lower limit on the 
plasma temperature kT c > 20 keV. Since the Kompaneet's 
equation (1) is valid only for kT c << m e c 2 , kT c = 25 keV 
has been chosen for the analysis (choosing kT c = 50 keV 
makes no significant qualitative difference in the results). 
These values of kT c and T give N esc = 1.44 (equation 2) 
and a Compton amplification factor A m 1.33. 

The model predictions and RXTE observations of the 
amplitude (root mean square, or RMS, fraction) as a func- 
tion of energy for the 830 Hz QPO (Berger et al. 1996) 
are compared in Figure 1, for different values of I and /. 



The overall normalization of the curve is fixed by choos- 
ing a best fit value of |<5T C |. The predicted amplitude is 
zero at E « 4 keV (comparable to the peak of the input 
spectrum) since the emergent spectrum pivots about this 
energy. This behavior is expected as oscillations in the 
soft photon Wien temperature and in the plasma temper- 
ature both preserve photon number. The lag is illustrated 
as a function of energy in Figure 2 (Vaughan et al. 1998) 
compared to the model prediction for the same values of I 
and / used in Figure 1. Reasonable fits to both variations 
are obtained for 4 km < I < 15 km and 0.5 < / < 0.85. 

Since the amplitudes of plasma and Wien temperature 
oscillations are nearly equal (i.e. |<5T C sa |<5T S |), the pho- 
ton fluctuations at high energies, Sn 7 H = Sn^(E « 15 
keV), are predominately caused by 5T C while the photon 
fluctuations at low energy, Sn^L = Sn 7 (E w 3 keV), are 
predominately caused by ST S . Thus, the time lag between 
5T C and Sti^h is Atu w Ni{t c , where Nh ~ 9 is the av- 
erage number of scatters a high energy photon undergoes. 
The time lag between ST C and Sn^L can be expressed as 
AtL ~ t c + At sc , where At sc is the time lag between (ST C ) 
and (ST S ). Thus the time lag between high and low energy 



photons is Ath s 



8t c 



At sc , so that when At sc > 8t c , 



the system will show soft lags. 

If the soft photon source is assumed to be a black body 
instead of a Wien peak, an additional ad hoc oscillation in 
the soft photon number flux is required for the predicted 
RMS fraction vs. energy to agree with observations. Thus, 
a Wien spectrum is favored for the soft photons, since it 
minimizes the number of oscillating components needed. 

In our model, the primary oscillation takes place in the 
plasma with the soft photon source responding to the vari- 
ation. Alternately, it is possible that the soft photon Wien 
temperature oscillates and, in order to keep the total lu- 
minosity of the plasma constant, induces variability in the 
plasma temperature. However, contrary to observations, 
such a model does not predict soft lags. Thus it seems 
that the QPO manifests itself only as an oscillation in the 
plasma temperature, rather than variations in other pa- 
rameters of the Comptonization model. 

The large value for the fraction of Comptonized pho- 
tons incident on the soft photon source, / ^ 0.5, suggests 
that the geometry of the hot plasma can be described as 
a corona on top of a cold accretion disk. The size of the 
corona (I « 5 km) indicates that it is a significant fraction 
of the radius of the disk. This result implies that the nearly 
coherent variability making up the kilohertz QPO is not 
limited to a narrow region and hints at the existence of a 
global mode in the corona. In such an interpretation, the 
high Q values (~ 100) characteristic of the high frequency 
oscillations would provide a strong constraint on the ex- 
citation and damping of the mode (Nowak et al. 1997). 
We point out that the Comptonizing model does not pro- 
vide insight into the the detailed geometry of the source. 
An alternative geometry would be an X-ray source (i.e., 
the neutron star) surrounded by a flaring accretion disk. 
This geometry introduces an additional time lag due to 
the light travel time between the Comptonizing plasma 
and soft photon source. 

Since the energy spectra did not exhibit an observ- 
able cutoff, the plasma temperature could be significantly 
larger than the assumed 25 keV. In that case, the aver- 
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age energy transfer per scatter would become a significant 
fraction of the photon energy and the Kompaneet's equa- 
tion, and the model predictions with it, would no longer 
be accurate (Katz 1987). Hence, simple extensions of our 
model to higher plasma temperatures are not useful, and 
detailed fitting to the observations will be warranted only 
when models incorporating non-standard geometries in so- 
phisticated time-dependent Monte Carlo methods (incor- 
porating energy balance) are developed. In this light, the 
ability of the simple model presented here to explain the 
main features of a kilohertz QPO's energy dependent prop- 
erties is highly encouraging. 



In the future, this work will be extended to other sys- 
tems exhibiting QPOs, particularly to the upper kilohertz 
QPO of the doublet observed in neutron star systems. 
Such an endeavor will not only check the robustness of the 
results presented here, but will also provide further con- 
straints on theoretical models of the QPO phenomenon. 
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Fig. 1. — Fractional RMS Amplitude vs. Energy. The model predictions for I = 5 km (solid curves) and I = 10 km (dashed curves) are 
compared to the data points from Berger et al. 1996. The vertical error bars denote one sigma errors, while the horizontal bars represent the 
energy bands used but do not account for the uncertainty in the PCA's energy resolution. Note that model curves 1 and 3 overlap. Curve 1: 
/ = 0.5, I = 5 kms and ST S = 0.585T c exp(i0.16ir). Curve 2: / = 0.75, I = 5 kms and 5T 3 = 1.75<5T c exp(i0.46ir). Curve 3: / = 0.5, I = 10 kms 
and ST a = 0A55T c exp(i0.27n). Curve 4: / = 0.75, I = 10 kms and ST a = 0.888T c exp(iOA8n). . 
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Fig. 2. — Relative Phase Lag (Radians) vs. Energy. The curves are model predictions for the same parameters as in figure 1. The data are 
taken from Vaughan et al 1998; as in Figure 1, the vertical bars are one sigma errors and the horizontal bars show the energy binning of the 
data. . 



